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The transient evolution of carriers in an intrinsic graphene under ultrafast excitation, which is 
caused by the collisionless interband transitions, is studied theoretically. The energy relaxation due 
to the quasielastic acoustic phonon scattering and the interband generation-recombination tran- 
sitions due to thermal radiation are analyzed. The distributions of carriers are obtained for the 
limiting cases when carrier-carrier scattering is negligible and when the intercarrier scattering im- 
poses the quasiequilibrium distribution. The transient optical response (differential reflectivity and 
transmissivity) on a probe radiation and transient photoconductivity (response on a weak dc field) 
appears to be strongly dependent on the relaxation and recombination dynamics of carriers. 



PACS numbers: 72.80.Vp, 78.67. Wj, 81.05.ue 

I. INTRODUCTION 

The transient response of photoexcited carriers un- 
der ultrafast interband pumping has been studied during 
the last decades in bulk semiconductors and heterostruc- 
tures (see Ref. 1 for review). The unusual transport of 
carriers in graphene is caused by a ncutrinolikc energy 
spectrum in gapless semiconductor, which is described 
by the Weyl- Wallace model 0], and a substantial mod- 
ification of scattering processes. Recently, the proper- 
ties of graphene after ultrafast interband excitation at- 
tract special attention. The experimental results in re- 
laxation dynamics ofphotoexcited electrons and holes 
were published in [!, U [f| and @ for epitaxial and exfo- 
liated graphene, respectively. The relaxation of nonequi- 
librium optical phonons, which are emitted by carriers af- 
ter photoexcitation, is studied in Ref. 7. The theoretical 
consideration of the carrier relaxation and generation- 
recombination processes caused by optical phonons is 
performed in @, Q . The quasielastic energy relaxation 
of carriers due to acoustic phonons was considered in 
[Tol . [Tl| for low energy carriers (at low temperatures or 
under mid-IR excitation). In particular, an interplay 
between energy relaxation and generation-recombination 
processes determines the relaxation dynamics of photoex- 
cited carrier distribution. [10[ To the best of our knowl- 
edge, both this interplay and the relaxation dynamics at 
low temperatures are not considered so far. Thus, the 
investigation of the transient response of carriers under 
these conditions is timely now. 

In this paper, we consider the transient response of 
an intrinsic graphene in case of ultrafast interband ex- 
citation in passive region, where the carrier energies are 
smaller than the optical phonon energy. Such a regime 
can be realized under the pumping in the mid-infrared 
(IR) spectral region or at low temperatures, when the 
peak of photoexcited carriers formed after the process 
of optical phonon emission, remains a narrow one. De- 



scribing the photoexcitation process, we restrict ourselves 
by the collisionless regime, when a pulse duration, r p , 
is shorter than the momentum relaxation time. Con- 
sidering the low-temperature transient dynamics of pho- 
toexcited carriers, one takes into account the intraband 
quasielastic energy relaxation due to acoustic phonons 
and generation-recombination interband transitions due 
to thermal radiation. The carrier-carrier scattering is 
described within two limiting regimes: (i) when the 
Coulomb interaction is unessential, and (ii) when inter- 
carrier scattering imposes the quasiequilibrium distribu- 
tion of carriers. With the obtained transient distribution 
of carriers, we analyze a time-dependent response on the 
probe field, i.e. we consider the transient reflection and 
transmission in THz and mid-IR spectral regions. The 
transient photoconductivity is also analyzed below, be- 
cause the energy relaxation corresponds to a nanosecond 
scale (the radiative recombination remains essential up 
to microsecond). 

Since the electron-hole energy spectrum and scattering 
processes are symmetric in an intrinsic graphene, the phe- 
nomena under consideration are described by the same 
distribution functions for electrons and holes, f p t- Such 
distribution is governed by the general kinetic equation 
El: 



^ = £-M/«|p} 



G{f\pt}, 



(1) 
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where the collision integrals Jk{ft\p} describe the re- 
laxation of carriers caused by the carrier-carrier scat- 
tering (k = cc), the acoustic phonons (k = ac), and 
the thermal radiation (k = r), respectively. The pho- 
togeneration rate, G{f\pt}, describes the interband ex- 
citation of electron-hole pairs by the mid-IR ultrafast 
pulse. Below Eq. (1) is solved with the initial condition 

fpt^-oc = fp eq ^ , where fp eq ^ is the equilibrium distri- 
bution. The transient response on a probe radiation is 
described by the dynamic conductivity due to interband 
transitions. The transient response on a weak dc field 
(photoconductivity) is considered with the use of the phe- 
nomenological model of momentum scattering suggested 
infill. 
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The analysis carried out below is organized as follows. 
The photoexcitation process under the inerband pump- 
ing is described in Sec. II. The transient evolution dis- 
tributions are given in Sec. Ill for the cases (?) and (ii). 
Section IV presents a set of results of transient reflectiv- 
ity and transmittivity, and also the transient photocon- 
ductivity. The discussion of the assumptions used and 
concluding remarks are given in the last section. Ap- 
pendix contains the microscopical evaluation of the in- 
terband photogeneration rate under ultrafast interband 
excitation. 



II. ULTRAFAST EXCITATION 

In the framework of the Weyl- Wallace model (spin- 
and valley-degenerate linear energy spectrum of carriers 
which is determined by the characteristic velocity %), 
the interband photoexcitation is caused by the in-plane 
electric field, W(Eexp(— i£lt) + c.c. where E is the field 
strength, u> is the frequency, and Wt is the envelope form- 
factor. Eq. (1) is transformed to the collisionless form 
on the initial intervals, when scattering mechanisms are 
not essential: df p t/dt = G{f\pt}. Using the boundary 
condition of Eq. (1), one can rewrite this equation in 
the integral form f pt = fj, eq) + f*^ dt'G{f\pt'}. The 
photogeneration rate here is evaluated in Appendix as 
follows: 



G{f\pt} = ( 6 ^ ) ir, j <It ic, : 



2v w p 



Q T 



(1 _ 2.fpt+T) ! 



(2) 



where the Pauli blocking factor (1 — 2f pt + T ) is responsible 
for the coherent Rabi oscillations of the excited carriers. 

Introducing the dimcnsionlcss intensity, I ex — 
(eEr p vw /hft) 2 , we consider below the linear regime of 
excitation which takes place if I ex <C 1 and f p t <C 1, so 
that the Pauli factor can be neglected (if Ml compara- 
ble to the equilibrium temperature T one has to use the 
equilibrium Pauli factor in Eq. (2)). Using the Gaussian 



form-factor w t 



2/7rexp —(t/r p ) with the pulse 



duration t p , [l4j one obtaines the photoexcited distribu- 



tion in the form 



Aex) 

Jpt 



I ex / dt'wt' I drwt'+r cos 



2vw5p 



(3) 



Here bp = p — pn is centred in the characteristic momen- 
tum = hfl/2vw- The evolution of photoexcited distri- 
bution, f p t /lex, is shown in Fig. 1. The distribution is 
dependent on t/r p and bp) Ap, where Ap = Ti/2v\yT p de- 
termines the width of distribution which is proportional 
to r~ x . For t t p , the integrations in Eq. (3) can be 



exactly performed and we obtain steady-state distribu- 
tion after the photoexcitation pulse, f p ex ^ = f^l^, as 
the Gaussian peak of width oc Ap: 



f(ex) 
Jp 



77 t ~(Sp/V2Ap) 2 
2 A ex c- 



(4) 



Thus, at t > 2t p (see Fig. 1) one can omit the photogen- 
eration rate in Eq. (1) using instead the initial condition: 



Ut= = 4 e9) + f P 



p(ex) 



(5) 



which is given as a sum of the equilibrium and photoex- 
cited contributions. The condition (5) can be used di- 
rectly in case of weak intercarrier scattering. In case of 
optical excitation, with a subsequent emission of cascade 
of 2Af optical phonons of energy hujo, the photoexcited 
distribution can be written in the form (5) where 5p is 
centred in pu = (HQ — 2J\fhio )/2vw and Ap is included 
an additional broadening during the cascade emission. 




t/r 



FIG. 1: Temporal evolution of photoexcited distribution f£f 
normalized to I ex versus dimensionless momentum and time, 
Sp/2Ap and t/r p . 

Under an effective intercarrier scattering, one needs 
to calculate the initial temperature and concentration of 
carriers. The photoexcited concentration and energy of 
carriers, which are described by the peak of distribution 
(4) are given by 



An ex 
AE„ X 



4 

1? 



p 



1 

v w p 



r 



(ex) _ Iex(uj/Vw) 2 



2u>T n 



1 

ftuj/2 
(6) 

where L 2 is the normalization area. One obtains 
AE ex /n ex = hZJ/2, for the Gaussian shape of pulse, i.e. 
the averaged energy per generated particle is equal to 
the excitation energy. In case of optical excitation, with 
J\f optical phonons emitted, the energy per photoexcited 
particle, AE ex / An ex , agrees closely with huj—JVhujQ (see 
above). 
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FIG. 2: Initial maximum distribution (a) and effective tem- 
perature (b), fex and T ex , versus pumping (An el /n T cc I ex 
for Tluj — 60 meV and 120 meV (solid and dashed curves, 
respectively) . 



If t p -C t cc <C T ac ,r , where t cc , T ac , and r r correspond 
to the intercarrier scattering, the energy relaxation, and 
the generation-recombination processes, respectively [the 
Coulomb-controlled case (ii)], the dominanting carrier- 
carrier scattering imposes the quasi-equilibrium distri- 
bution 



I pi 



( v w p ~ IH . , . 
exp \ — | + ! 



(7) 



III. ENERGY RELAXATION AND 
RECOMBINATION 

In this section we analyze the transient evolution of 
f p t caused by the energy relaxation and recombination 
processes. We consider the cases (i) and (m), when the 
initial condition is given by Eq. (5) and written through 
/ex and T ex plotted in Fig. 2. 



A. Weak intercarrier scattering 

If the carrier-carrier scattering is ineffective [case (i)], 
the distribution f pt is governed by the kinetic equation 
(1) without the cc-contribution 



dfpt 

at 



S(qe) 



p 2 dp 



dfpt 
dp 



/ P t(l — fpt) 



PT 

2 fpt) 



(9) 



and with the initial condition (5) used instead of gen- 
eration rate. Here we substituted the explicit expres- 
sions of the collision integrals for the quasielastic acoustic 
scattering approximation (written in the Fokkcr-Planck 
form) and for the generation-recombination processes, 
see discussion in [l(| ■ The Planck distribution N 2p /p T 
is written through px = T/vw while the energy relax- 
ation rate v v qe ^ = v qe p/% and the rate of radiative transi- 
ts) 

tions Vp — v r p/h are written through the characteristic 
velocities v qe oc T and v r [l5j . 

The boundary conditions are imposed by both the con- 
dition /p_+oot = 0, which is transformed into the require- 
ment 



with the effective temperature T t and the quasichemi- 
cal potential /i t . If t cc <C t <C T afvr , the initial values 
T ex — T t —>o and f ex = / p =ot->o are determined from the 
concentration and energy conservation requirements: 



2 00 

- ( -r^- I / dxxf x 
7r \hvwj J 



1 



TlT 

Ex 



l\n ex 



(8) 



where the function f x is introduced according to f x = 
fex/[e x (l ~ fex) + fex}- Using An ex and AE ex given by 
Eq. (6) and solving the transcendental system (8) one ob- 
tains the initial values f ex and T ex . The calculations here 
and below are performed for the nitrogen temperature, 
T =77 K, the excitation energies 2vwPu =120 meV (CO2 
laser) and 60 meV (as an example of interband excita- 
tion with subsequent optical phonon emission), and the 
broadening energy h/r p ~6.6 meV, which corresponds to 
the pulse duration ~0.1 ps. In Fig. 2 we plot f ex and 
T ex versus the pumping level which is proportional to 
An ex /nT- Fast increase of T ex and fast decrease of f ex 
take place for An ex /riT < 1, while a linear increase of 
these values are realized if An ex /nT > !• 



df 



pi 



dp 



fpt 

PT 



< const, 



(10) 



and Eq. (9) at p = which is transformed into the ini- 
tial condition / p= ot = 1/2 + fp=o exp[—(y r /vw)Tt/K]. 
According to Eq. (4) one obtains /^ffg = 
^Jir/2I e x exp[— (f2r p ) 2 /2] <C 1 and one can neglect the 
second contribution in this initial condition, so that 
f P =ot — 1/2- Numerical solution of the Cauchy prob- 
lem given by Eqs. (5), (9), and (10) is obtained below by 
the use of the iteration procedure. [l6| 

In Fig. 3 we demonstrate the evolution of the distri- 
bution fpt at 77 K for the cases when carriers are ex- 
cited around the energies 60 meV and 30 meV. The delay 
times are marked in panels a-c and the pumping levels 
are determined through the initial peak value, given by 

it I II ex, see Eq. (4)). Under mid-IR pumping with 
pulse duration t p =0.1 ps and the spot sizes ~0.5 mm 
the above-used pumping levels correspond to the pulse 
energies ^85 pJ and ~17 pJ for Figs. 3a and 3b, re- 
spectively, see 17] for experimental details. Under op- 
tical pumping (htt ~1.6 eV) and subsequent emission 
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take forms: [16 




20 40 60 

v w p (meV) 

FIG. 3: Distribution f p t versus carrier energy pvw for dif- 
ferent delay times (marked) and excitation conditions: (a) 
hZJ =120 meV and I ex =0.26, (b) hZJ =120 meV and 
lex =0.05, and (c) Tuj =60 meV and I ex =0.21. 



of phonon cascade, the pumping level in Fig. 3c corre- 
sponds to the pulse energy ~12 nJ (duration and size are 
the same as above). One can see that the transient evo- 
lution of distribution occurs in two stages: energy relax- 
ation and recombination. During the first stage (about 
t <50 ns, which is dependent on position and maximum 

value fp CX ^] compare with Figs. 3a-c) the initial peak 
is tranformed into the quasiequilibrium high-energy tail 
(with the equilibrium temperature caused by the energy 
relaxation) which is connected to the low-energy equilib- 
rium distrbution. During the next stage (up to 1 /is) the 
high-energy tail shifts to the lower energies and trans- 
forms into the equilibrium distribution due to effective 
radiative recombination in low-energy region. 



B. Coulomb-controlled case 

In the carrier-carrier scattering case (ii), one has to 
describe the transient evolution of the effective temper- 
ature T t and the maximum distribution f t = f p =ot, 
that replaces the chemical potential. Since the in- 
tercarrier scattering change neither the concentration, 
n t = (4/L 2 ) X] p nor the energy of carriers, E t = 
(4/ L 2 )J2 p v wPfpt, the balance equations for n t and E t 



d_ 

dt 



m 



L 2 ^ 
p 



MMp} 

VWP [Jac{ft\p} + Jr{ft\p}] 



(11) 



Further, we transform the balance equations, expressing 
the left-hand side of (11) through T t and f t as follows: 



dt 



R 



(i) 



(12) 



dt 



(tM 2) ) 



R[ 2) + Qt- 



Here the coefficients A 



(1,2) 



are written as A 



(?) 



/ °° dxx q f x t, where the quasiequilibrium distribution is 

given by f xt = f t / [e x (l - f t ) + f t ], so that A ( t q) /T l t are 
only depend on f t . After substitution of the collision 
integrals J r [Tol . [l8| and integration, the generation- 
recombination contributions to Eq. (12) are obtained 
in the form 



R 



(9) 



2v r T t q+ 



dxx q+ f xt 



e 2x (l~ft) 2 

e x2T t /T _^f2 



- 1 



(13) 

Similarly, the energy relaxation contribution is written 
by the use of J ac as follows 



Qi 



T-T t v qe T t A 
T vwh 



dxx e f xt 



(14) 



The initial conditions for the system (11) are written as 
Tt=o = T ex and ft=o — fex- 

Numerical solution of the nonlinear system (11) is per- 
formed using the iteration procedure. In Fig. 4 we plot 
Tt/T and ft versus time. Temperature relaxes to the 
equilibrium one during the energy relaxation times (^ 
100 ns) while ft, which is determined by the chemical 
potential fit, relaxes to 1/2 over 1 [is (the recombination 
time scale), in analogy with the case (i). Notice, that af- 
ter the fast energy relaxation, one obtains ft >l/2 [dot- 
ted line in Fig. (4b)], i.e. the low-energy electron-hole 
pairs appear to be unstable. [l9[ 

Fig. 5 shows the plot of temporal evolutions of the 
energy per particle and concentration, E t / n t and nt [see 
the definitions before Eq. (11)], for the cases (i) and (ii). 
The relaxation processes to the equilibrium (at nitrogen 
temperature, £ , f^ 00 /nt^ o —14.5 meV and nt-+oo — 5.3- 
10 9 cm -2 ) occur during the same scales as in Figs. 3 and 
4. The temporal dependencies of nt obtained for both 
cases are in good agreement (the carreir-carrier scattering 
does not change concentration) while E t demonstrates 
a different evolution for cases (i) and (ii) at t <50 ns. 
This is because of drift and decrease of photoexcited peak 
during the energy relaxation time, see Fig. 3. 
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FIG. 4: Temporal evolution of effective temperature, Tt (a), 
and maximum distribution, ft (b), for different excitation con- 
ditions: (1) hZJ =60 meV and I ex =0.21, (2) hZJ =60 meV and 
lex =0.1, (3) fu3 =120 meV and I ex =0.052, and (4) Ku =120 
meV and I ex =0.026. 




0.4 0.6 

f Qis) 

FIG. 5: Energy per carrier (a) and concentration (b) versus 
time. Solid and dotted curves correspond to the cases (i) and 
(ii) , respectively; excitation conditions (1) - (4) are the same 
as in Fig. 4. 



IV. TRANSIENT RESPONSE 

Here we turn to consideration of the response of pho- 
toexcited carriers on a probe radiation (reflection and 
transmission in the THz and mid-IR spectral regions) 
and on a weak dc electric field (photoconductivity). The 



transient electrodynamics of graphene is described using 
the time-dependent dynamic conductivity, <7 wt , which is 
caused by the collisionless interband transitions, see Ap- 
pendix B. The transient photoconductivity is calculated 
by the use of the phenomenological model of momentum 
relaxation suggested in fljjj ]. 



A. Reflection and transmission 

To calculate the transient reflectance and transmit- 
tance of the graphene sheet placed at z — on the 
in-plane electric field E zt exp(— iwt) propagated along 
0Z, we apply the wave equation, see [20j and references 
therein. The induced current density, cr wt _E z= o, is lo- 
cated around z = and direction of in-plane field E zt 
is not essential due to the in-plane isotropy of the prob- 
lem. Separating the incident radiation, E in e " z , with 
the wave vector k w = ui/c, we write the field distribution 
outside of the graphene sheet in the form: 



E z t — 



E», 



E 



Z < 
Z > 



(15) 



where — •Jeui/c is the wave vector in the substrate 
with the dielectric permittivity e. The transmitted and 



reflected electric fields, and E\' 1 , are determined 

from the boundary conditions at z — > as follows: 



(r) 



El 



(*) 



E, 



(t) 



E,, 



1 + 4, 



E-,, 



1-4, 
1 + 4 



(16) 



Here we introduce the dimensionless factor A w t — \fe + 
(47r/c)<7 wt . The reflection and transmission coefficients, 



= \E[ r) \ 2 
through A u t according to 

i — \ ,, 

Rut = 



FU = \En 2 /Ef n and T L 



= \E x t L '\' z /Ef n , are written 



1 + 4 



Tut — 



U + 4 



(17) 



Using a ut determined by Eqs. (B3) and (B4), we con- 
sider below the differential changes in reflectivity and 
transmissivity, (AR/R) u t = {R^t — ) / and 
(AT/TV = {T ut - Ti eq) )/Ti eq \ which are written 
through the equilibrium reflection and transmission co- 
efficients, R ( J q) and T^ eq) . 

The evolution of the differential reflectivity for the 
cases (i) and (ii) are shown in Figs. 6a and 6b, 6c, 
respectively. If the Coulomb scattering is not effective 
[case (i)], the distribution of carriers relaxes during the 
energy relaxation time scale (around 10 ns, cf. with Fig. 
3), when a quenching of photoexcited peak takes place (if 
hcu is comparable with the peak energy). In case (ii) any 
peculiarities of the spectrsal dependencies at stort times 
are absent because the initial distribution is transformed 
into the quasiequilibrium one during times ~ r cc — > 0. 
The further evolution of (AR/R)^ is limited by the 
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FIG. 6: (a) Spectral dependencies of differential reflectivity, 
(AR/Kjut, for different delays (marked) at the excitation con- 
ditions: (a) hZJ =120 meV and I ex = 0.052 in the case (i), 
(b) hU = 120 meV and I ex =0.052 in the case (ii), and (c) 
hZJ =60 meV and I ex =0.21 in the case (ii). 



generation-recombination process and extended up to mi- 
croseconds. In the THz spectral region (Tllo >10 meV is 
considered here because we neglect the intraband relax- 
ation), the differential reflectivity increases and changes 
a sign. In the high-energy region, (AR/R) ut decreases 
monotonically with oj and t and does not exceed ~ 
for the near-IR spectral region. Beside of this, the re- 
sponse is approximately proportional to the pumping in- 
tensity, I ex , and (Ai?/i?) wt increases with increasing of 
the photoexcitation energy, TiUj (cf. Figs. 6b and 6c). 



In Fig. 7 we plot the differential transmissivity for the 
cases (i) and (ii) under the same excitation conditions. 
Once again, in the high-frequency region the differen- 
tial transmissivity decreases slowly (during a microsec- 
ond time scale) and (AT/T% t does not exceed ~ 10~ 4 
for the near-IR spectral region. In the THz spectral re- 
gion, (AT/T) wt increses and changes the sing in the same 
manner as (AR/R) ut (cf- Figs. 6 and 7). The depen- 
dencies on the excitation parameters (I ex and Tuj) are 
also similar to the reflectivity. Additionally, in case (i) 
a fast (at t <10 ns) quenching of the photoexcited peak 
contribution in the spectral region ~ hfl takes place. 
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FIG. 7: (a) Differential transmissivity, (AT/T) u t, versus fiu) 
and t for cases (i) and (ii) [panels (a) and (b), respectively] 
at the same excitation conditions: Tiuj =120 meV, I ex =0.052. 



B. Photoconductivity 

Finally, we consider the transient photoconductivity, 
i.e. the response of the photoexcited carriers to the weak 
dc electric field. Since the momentum relaxation is gov- 
erned by elastic scattering mechanisms, [HI one can use 
the following expression for the dc conductivity at- 



ft = fo 



Ir 



dp/-. 



pi 



(18) 



Here l c is the correlation length characterizing the disor- 
der scattering and the function *&(z) = e~ z Ii(z 2 )/z 2 is 
written through the first order Bessel function of imagi- 
nary argument, I\(z). The normalized conductivity, <To, 
is introduced for the case of short-range scattering, when 
l c = 0. The distribution / p= ot is shown in Fig. 4b for 
the case (ii) while 2/ p=0 t = 1 for the case (i). If l c = 0, 
one obtains at/era = 1, i.e. there is no transient photo- 
conductivity for the case (i); for the case (ii) one obtains 
ft/fo = 2/t and the transient photoconductivity is clear 
from Fig. 4b. 

If l c ^ 0, the transient evolution of conductivity is 
shown in Fig. 8. For the definiteness, it was assumed 
that Z c =10 and 30 nm and variations of <j t are increased 
with l c essentially due to contribution of high-energy car- 
riers. Similar to Sec. IVA, one can separate two stages 
of evolution: the fast decrease of ct* due to energy relax- 
ation (up to ~ 30 4- 50 ns for the conditions considered) 
and the slow quenching of a t due to carrier recombina- 
tion. If t > 1 /is, the conductivity approaches to the 
equilibrium values: at^oo/&o =1.445 if Z c =30 nm and 
ft— >oo/fo =1-035 if Z c =10 nm. Since the transient con- 
ductivity can be measured for the subnanosecond time 
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3 2 4 I =30 nm 




FIG. 8: Temporal evolution of conductivity for excitation con- 
ditions (l)-(4) which are the same as in Fig. 4 for the corre- 
lation length l c =30 nm and 10 nm. Solid and dashed curves 
are correspondent to the cases (ii) and (i). 



scale [21( , such a scheme can be used for verification both 
energy relaxation and recombination mechanisms. 



V. CONCLUDING REMARKS 

To summarize, we have considered both the interband 
ultrafast photoexcitation and the relaxation dynamics of 
the carriers in an intrinsic graphene. In contrast to the 
measurements [1, 13, S @ an d calculations @, [§] per- 
formed, where the evolution corresponds to the subpi- 
cosecond time scales due to the opticlal phonon contri- 
bution, here we consider the slow relaxation of the low- 
energy carriers. The distribution of carriers at T =77 
K is obtained for the limiting cases with negligible or 
dominating intercarrier scattering when the energy relax- 
ation and generation-recombination processes are caused 
by the quasielastic acoustic phonon scattering and ther- 
mal radiation, respectively. The initial distribution is 
obtained in the framework of the linear, with respect to 
pumping, approximation for the collisionlcss regime of 
the interband transitions. The transient optical response 
on the probe radiation (transmission and reflection) as 
well as on the weak dc field (transient photoconductiv- 
ity) appears to be strongly dependent on the relaxation 
and recombination dynamics of carriers. 

Next, we discuss the assumptions made. The main re- 
strictions of the results presented are the consideration 
of the low-energy carriers, when the interaction with op- 
tical phonons is unessential, and the single generation- 
recombination mechanism (due to thermal radiation) is 
taken into account. These conditions are realized at low 
temperatures under the mid-IR ultrafast excitation [ItJ 
of the clean sample (e.g. suspended graphene [HI). Such 



an approach can be used for the case of optical inter- 
band excitation, when the low-energy initial distrbution, 
with a phcnomenological broadening, is formed after the 
cascade process of optical phonon emission. The con- 
sideration is restricted by the radiative recombination 
(the Auger processes are forbidden due to the symme- 
try of electron-hole states [HI), with the characteristic 
time scales up to microseconds. Any visible contribu- 
tion of other generation-recombination mechanism (e.g., 
because of disorder-induced interband transitions with 
acoustic phonons, or under intercarrier scattering) leads 
to fast decrease of photoresponse. Such a regime requires 
an additional investigation but the quasielastic energy re- 
laxation stage is described by the presented results. 

The rest of assumptions are rather standard. The con- 
sideration in Sec. Ill is limited by the simple cases (i) 
and (ii), with and without the intercarrier scattering. 
The main peculiarities of the response under considera- 
tion are similar for both cases but the complete descrip- 
tion of the nonequilibrium carriers had been performed 
neither under optical excitation, nor under high dc field, 
see [l8|, [24[ and Refs. therein. The description of the 
momentum relaxation in Sec. IV is based on the phe- 
nomenological model of Ref. 13. The utilization of the 
quasielastic energy scattering and the collisionlcss inter- 
band photoexcitation appear to be rather natural. The 
listed assumptions do not change either the character of 
the response or the numerical estimates. 

In closing, the peculiarities of the transient optical re- 
sponse (transmission and reflection) as well as of the tran- 
sient photoconductivity appear to be useful tool in or- 
der to verify the relaxation and generation-recombination 
mechanisms of carriers. Thus, in addition to the recently 
obtained experimental results 0, 0> IE @> measurements 
under mid-IR excitation and at low-temperature will be 
useful for characterization of graphene. 



APPENDIX A: GENERATION RATE 

Below we describe the interband carrier excitation un- 
der ultrafast mid-IR pumping ~E t exp(— iQt) for the col- 
lisionless case, when r p is shorter than relaxation times. 
The photogeneration rate into the a-state is based on the 
general expression (see [l| and Sec. 54 in Ref. 12) 



G a t = 2i?e 



(—)' 



a U hr ' h [(E t+T • v) , p t+T ] e- l ' hT ' n , (E t • v)~ 



a), 



where pt is the density matrix, v is the velocity operator, 
and A — * +0. Since the collisionlcss regime of photoexci- 
tation, we calculate (Al) with the use of the free states 
|Zp) and the energy e\ v where I = ±1 stands for c- or v- 
bands and p is the 2D momentum. Neglecting the non- 
diagonal components of the density matrix p t and using 
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the distribution functions /; p t, one obtains the genera- 
tion rate 

o 

G{f\lpt}= (^) 2 / dTe XT - inT e l ^- e -^ T / % 

— OO 

x(lp|(E t+r • v)| - lp)(-lp|(E t • v)+|lp) (A2) 

X (/-lpt+r - flpt+r) + C.C. , 

moreover G{f\ — Ipt} = —G{f\lpt} according to the 
particle concervation law. Next, we separate the envelope 
form-factor wt using E 4 = E,w t and take into account 
the in-plane isotropy of the problem, when one arrives 

to the averaged matrix element |(+lp |(E ■ v)| — lp) | 2 = 
(Evw) 2 /2. As a result, we obtain the in-plane isotropic 
generation rate G{f\pt} — ±G{/| ± Ipt} in the following 
form: 

— OO 

xel (2v wP )r/H (/ _ lpf+r _ /lpt+r) + c . c . . (A3) 

Finally, using the electron-hole representation and re- 
placing the filling factor here by (1 — 2f pt ), we arrive 
to Eq. (2). 

APPENDIX B: DYNAMIC CONDUCTIVITY 

The response of graphene on the in-plane probe field 
Eexp(— iut) is described by the dynamic conductivity 

turn 

^-^Ed- 2 W (bi) 

CU-U — 
P 

x ( 1 1 ^ 



with A — > +0. The parametric time dependency of <r ut is 
valid if the time scales under consideration exceed uj^ 1 . 
It is convenient to separate the time-independent con- 
tribution, a u , described the undoped graphene in the 
absence of photoexcitation, when f pt vanishes. Using 
the energy conservation law one obtains Rea^ = e 2 /47i. 
In the framework of the Weyl- Wallace model, the Im- 
contribution into W w appears to be divergent at p — > oo. 
It is convenient to approximate Ima 7 ^ as a sum of oc uj^ 1 
and oc u) terms, which correspond to the contributions of 
the virtual interband transitions and the ion background, 
correspondingly. As a result, we obtain: 

where the characteristic energies, e m ~0.1 eV, and ~ 
6.8 eV are correspondent to the recent measurements of 
the graphene optical spectrum. [26] 

Next, substituting the time-dependent distribution f pt 
obtained in Sec. Ill into the dynamic conductivity (Bl) 
one transforms the real and imagional parts of a^t as 
follows 

e 2 

Rea wt = — [l-2F{ Pui ,t)], (B3) 

OO 

_ e 2 f dyy 2 

haa ut = Imer w -V / ?F(p u y,t). 

nh J l-y 2 
o 

Here V means the principal value of integral. We also 
introduced the function F(p, t) = f pt for the case (i) and 

F fay, t) = /; — — (B4) 

exp[(huj/Tt)y][l ~ ft) + ft 

for the case (ii), when a^t is determined both the effective 
temperature and the carrier concentration, T t and ft- 



[1] J. Shah, Ultrafast Spectroscopy of Semiconductors and 
Semiconductor Nano structures, (Springer, New York, 
1996); F. T. Vasko and A. V. Kuznetsov, Electron States 
and Optical Transitions in Semiconductor Heterostruc- 
tures (Springer, New York, 1998). 

[2] E. M. Lifshitz, L. P. Pitaevskii, and V. B. Berestetskii, 
Quantum Electrodynamics, (Butterworth-Heinemann, 
Oxford 1982); P. R. Wallace, Phys. Rev. 71, 622 (1947). 

[3] J. M. Dawlaty, S. Shivaraman, M. Chandrashekhar, F. 
Rana, and M. G. Spencer Appl. Phys. Lett. 92, 042116 
(2008). 

[4] D. Sun, Z.-K. Wu, C. Divin, X. Li, C. Berger, W. A. de 
Heer, P. N. First, and T. B. Norris, Phys. Rev. Lett. 101, 
157402 (2008). 

[5] P. A. George, J. Strait, J. Dawlaty, S. Shivaraman, Mvs. 
Chandrashekhar, F. Rana, and M. G. Spencer, Nanolet- 



ters 8, 4248 (2008). 
[6] R. W. Newson, J. Dean, B. Schmidt, and H. M. van Driel, 

Opt. Exp. 17, 2326 (2009). 
[7] H. Wang, J. H. Strait, P. A. George, S. Shivaraman, V. B. 

Shields, Mvs Chandrashekhar, J. Hwang, F. Rana, M. G. 

Spencer, C. S. Ruiz- Vargas, and J. Park. larXiv:0909.4912l 
[8] S. Butscher, F. Milde, M. Hirtschulz, E. Malic, and A. 

Knorr, Appl. Phys. Lett. 91, 203103 (2007). 
[9] F. Rana, P. A. George, J. H. Strait, J. Dawlaty, S. Shiv- 
araman, Mvs Chandrashekhar, and M. G. Spencer, Phys. 

Rev. B 79, 115447 (2009. 
[10] F. T. Vasko and V. Ryzhii, Phys. Rev. B 77, 195433 

(2008); A. Satou, F. T. Vasko and V. Ryzhii, Phys. Rev. 

B 78, 115431 (2008). 
[11] R. Bistritzer and A. H. MacDonald, Phys. Rev. Lett. 

102, 206410 (2009). 



9 



[12] F. T. Vasko and O. E. Raichev, Quantum Kinetic Theory 
and Applications (Springer, New York, 2005). 

[13] F. T. Vasko and V. Ryzhii, Phys. Rev. B 76, 233404 
(2007). 

[14] The form-factor w t is normalized according to the con- 
dition dtuit = T p . The shape of normalized form- 
factor has little effect on the transient photoexcitation 
under consideration because the ultrafast response is de- 
termined fundamentally by the pulse duration, r p . 

[15] Here we use the characteristic velocities v ac ~2.5xl0 
cm/s (for the nitrogen temperature) and v r —41.6 cm/s 
(for graphene sheet placed between Si02 substrate and 
cover layer), see [To( ] . 

[16] D. Potter, Computational Physics (J. Wiley, London, 
1973). 

[17] T. Elsaesser and M. Woerner, Physics Reports 321, 253 
(1999). 

[18] O. G. Balev, F. T. Vasko and V. Ryzhii, Phys. Rev. B 

79, 165432 (2009). 
[19] According to Eq.(B3), the negative interband absorption 

(oc Recr^t) takes place, if f pt > 1/2, see [T^] and Refs. 

therein. This low-energy instability is suppressed due to 

an effective intraband (Drude) absorption. 
[20] L. A. Falkovsky, Phys. Usp. 51, 887 (2008); T. Stauber, 



N. M. R. Peres, and A. K. Geim, Phys. Rev. B78, 085432 

(2008) ; M. V. Strikha and F. T. Vasko, submitted. 

[21] T. Yao, K. Inagaki, and S. Maekawa, in Proceedings of the 
11th International Conference on the Physics of Semicon- 
ductors (Polish Scientific Publishers, Warszawa, 1972), 
Vol. 1, p. 417. 

[22] G. Li, A. Luican, and E. Y. Andrei, Phys. Rev. Lett. 102, 
176804 (2009); P. Neugebauer, M. Orlita, C. Faugeras, A. 
L. Barra, and M. Potemski, Phys. Rev. Lett. 103, 136403 

(2009) . 

[23] M. S. Foster and I. L. Aleiner, Phys. Rev. B 79, 085415 
(2009). 

[24] A. Akturkaand N. Goldsman, J. Appl. Phys. 103, 053702 

(2008) ; R. S. Shishir and D. K. Ferry, J. Phys.: Condens. 
Matter, 21, 344201 (2009). 

[25] R. R. Nair, P. Blake, A. N. Grigorenko, K. S. Novoselov, 
T. J. Booth, T. Stauber, N. M. R. Peres, and A. K. Geim, 
Science 320, 1308 (2008); T. Stauber, N. M. R. Peres, 
and A. K. Geim, Phys. Rev. B78, 085432 (2008); K. F. 
Mak, M. Y. Sfeir, Y. Wu, C. H. Lui, J. A. Misewich, and 
T. F. Heinz, Phys. Rev. Lett. 101, 196405 (2008). 

[26] M. Bruna and S. Bonilla, Appl. Phys. Lett. 94, 031901 

(2009) . 



